arXiv:1507.05237v2 [astro-ph.GA] 16 Dec 2015 


Draft version December 17, 2015 

Preprint typeset using DT^X style emulateapj v. 01/23/15 


SENSITIVE 21CM OBSERVATIONS OF NEUTRAL HYDROGEN IN THE LOCAL GROUP NEAR M31 

Spencer A. Wolfe 

Dept, of Physics & Astronomy, West Virginia University, Morgantown, WV 26506 


Felix J. Lockman 

National Radio Astronomy Observatory 1 , Green Bank, WV 24944 


D.J. Pisano 

Dept, of Physics & Astronomy, West Virginia University, Morgantown, WV 26506 
Draft version December 17, 2015 

ABSTRACT 

Very sensitive 21cm HI measurements have been made at several locations around the Local Group 
galaxy M31 using the Green Bank Telescope (GBT) at an angular resolution of 9fl, with a 5cr detection 
level of Nhi = 3.9 x 10 17 cm -2 for a 30 km s -1 line. Most of the HI in a 12 square degree area almost 
equidistant between M31 and M33 is contained in nine discrete clouds that have a typical size of a 
few kpc and HI mass of 10 5 M 0 . Their velocities in the Local Group Standard of Rest lie between 
— 100 and +40 km s -1 , comparable to the systemic velocities of M31 and M33. The clouds appear 
to be isolated kinematically and spatially from each other. The total HI mass of all nine clouds is 
1.4 x 10 6 M 0 for an adopted distance of 800 kpc with perhaps another 0.2 x 10 6 M 0 in smaller clouds 
or more diffuse emission. The HI mass of each cloud is typically three orders of magnitude less than 
the dynamical (virial) mass needed to bind the cloud gravitationally. Although they have the size and 
HI mass of dwarf galaxies, the clouds are unlikely to be part of the satellite system of the Local Group 
as they lack stars. To the north of M31, sensitive HI measurements on a coarse grid find emission that 
may be associated with an extension of the M31 high-velocity cloud population to projected distances 
of ~ 100 kpc. An extension of the M31 high-velocity cloud population at a similar distance to the 
south-east, toward M33, is not observed. 

Keywords: galaxies: halos - galaxies: ISM - intergalactic medium - Local Group 


1. INTRODUCTION 

The two largest galaxies in the Local Group, M31 and 
the Milky Way, have a substantial amount of gas residing 
in a circumgalactic medium (CGM, also called a gaseous 
halo), outside of their disks. Their CGM is dominated 
by ionized gas, but also contains neutral hig h-veloci ty 
clouds (HVCs) observed in the 21cm line 


Sembach et al.|2Q 03; Thil ker et al.1 2004; West meier et al 


2008 Shull etal. 


I2009H . 

,s likeh 


Putman et al.|2009 Lehner et al 


2012| |2015|). Gas likely associated with M31 is seen m 

absorption against backgrou nd AGN to project ed dis¬ 
tances of at least 300 kpc (Lehner et al. 2015). The 
Milky Way may have a similar CGM less easily separated 
from disk gas because of projection effects, but manifest 
in the stripping o f gas from dw arf sp he roidals at dis- 

Spekkens 
M31 does 

extend this far, it encompasses the smaller spiral M33, 
which itself has a modest population o f neutral HVCs 


tances to 300 kpc (IGrcevich fe Putman 
let al.||2014[ |Gatto et al. 112013 


2009 
11 the CGM of 


(Grossi et al. 2008 Putman et al. 2009). The CGM of 


the Milky Way contains the Magellanic Stream (MS) 
which extends at least 200° across the sky in HI emis¬ 
sion and whose mass is probably dominated by ionized 
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1 The National Radio Astronomy Observatory is a facility of the 
National Science Foundation operated under a cooperative agree¬ 
ment by Associated Universities, Inc. 


gas (|Nidever et al.||2010| |Fox et al.||2014 ). 

There are also neutral atomic hydrogen (HI ) clouds 
in the Local Group whose connection with individual 
galaxies is not understood. Compact high velocity 
clouds (CHVCs) and ultra-compact high velocity clouds 
(UCHVCs) are of small angular size and relatively iso¬ 
lated, and are candidates for low mass galaxies that may 


lack star fo rmation entirely (de Heij et al. 


et al. 


et al. 


2013). For a variety of reasons, CL 


2002 


Adams 


VUs are now 


thou ghtt o reside in the Milky Way CGM (e.g., Sternberg 


2002), but the location and nature of the recently^ 


discovered UCHVCs is less clear. 

In the ir study of 21cm emiss ion from HI in the Local 
Group, Braun & Thilker ([2004 hereafter BT04), discov¬ 
ering 


ered extended regions of H1 around the galaxy M31 that 
formed a partial bridge to the galaxy M33. This emission 
was detected at the extremely low levels of Nhi^ 10 17 
cm -2 , about two orders of magnitude below the typical 
column density detectab le in extragalactic 21cm obser¬ 
vations (Heald et al.| 2011). The BT04 observations were 
made with a rather coarse angular resolution of 49' and 
the origin of the neutral gas was uncertain, but the dif¬ 
fuse HI appeared to connec t the system ic heliocentric 
velocities of M31 and M33 ((Lewis et al.|[20131). BT04 
proposed that the HI arose from condensation in a dark 
matter-dominated filament connecting the two galaxies. 
Another suggestion was that i t resu l ted fr om a tidal en¬ 
counter between the galaxies (Bekki 2008). Subsequent 
observations of part of the region using the Green Bank 
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Telescope (GBT) at 9' angular resolution confirmed the 
reality of the emission, though those data lacked thesen- 
sitivi ty to reveal any detailed structure (Lockman et al. 
2012). As the existence and properties of CGM and 
intra-group gas is critical to our understanding ofthe for¬ 


mation and evolution of galaxies (e.g. Dave et al. 2001 


_ galaxies (e.g. _ 

Fnkngita fc Peebles|2006l|Chen et al.|20lO[ Butman et al 


2012[|Gen|2013||Leh ner et al. 2015), we have undertaken 

a major survey of the area around and between M31 and 
M33 using the Green Bank Telescope, which provides 
both the sensitivity needed to detect this extremely faint 
emission, and the angular resolution to discern some of 

its structure. _ 

In a previous paper (Wolfe et al.|2013 hereafter Paper 
I) we presented results of the first part of the study, which 
showed that a significant fraction of the HI detected in 
a 12 square-degree field south-east of M31 in the direc¬ 
tion of M33 arose in discrete structures which, assuming 
they are 800 kpc distant, have the size of dwarf galax¬ 
ies but no de tecta ble stellar component. Subsequently, 
Martin et al. (2013) suggested that there may be a stellar 
overdensity in the direction of one of the clouds, but the 
association between the gas and stars is not established 
at this time. The clouds have velocities similar to the 
systemic velocities of M31 and M33 and are thus prob¬ 
ably not part of the high-velocity cloud (HVC) system 
of either galaxy. Here we report on additional observa¬ 
tions that reveal more clearly the structure and content 
of the HI clouds, as well as new observations of selected 
directions to the north of M31. The distances to M31 
and M33 are ^ 750 and ^ 850 kpc, respectively (e.g. 


2012 Bono et al. 2010) so we assume that 
located between them lies at a distance of 


Riess et al 
the materia" 

800 kpc. 

2. OBSERVATIONS 

All data used here were obtained with the 100-meter 
Robert C. Byrd Green Bank Telescope (GBT, Prestage 


et al. 2009) with the dual-polarization L-bana receiver 
that has a total system temperature of < 18 K at ele¬ 
vations > 20°. At the frequency of the 21cm line, the 
telescope has a half-power beam-width of 9.'1. Spectra 
were measured using the GBT Spectrometer, which pro¬ 
vides a total velocity coverage > 1000 km s -1 in the 
21cm line at a channel spacing of 0.32 km s -1 . Spec¬ 
tra were calibrated using obser vations of 3C48 and th e 
antenna response analysis from Boothroyd et al. (2011). 
Velocities were measured with respect to the Local Stan- 
dard of Rest (LSR), but in this part of the sky, Vlsr 
differs from heliocentric velocities by only a few km s -1 . 

Figure [l] shows the observed areas on the BT04 
HI map oF the region. Data reduction was done with 
speci al procedures written in GBTIDL (Marganian et al. 
2006), the stray radiation corre ction and calibration fol¬ 
lowed Boothroyd et al. (2011), and special procedures 


were written to remove instrumental baselines from spec¬ 
tra in the data cube. 


2.1. The M31-M33 Field 

As described in Paper I, we used the GBT to map a 
6° x 2° field in right ascension and declination centered at 
J2000 = 01 h 16 m + 37°00 / . The observations were made 
while moving the telescope in right ascension at a fixed 
declination, binning the data in 3-second samples every 


17 18 19 20 21 
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Figure 1 . Locations of the current observations marked on a map 
of the total Nhi from the BT04 survey. The galaxies M33 and 
M31 lie at the lower left and upper right, respectively. Contours 
are for log(N m ) = 17.0, 17.3, 17.7, 18.0, 18.3, 18.7, 19.0, 19.3, 
19.7, 20.0, 20.3 and 20.7 (cm -2 ). Our “M31-M33” field, outlined 
with a box, lies between the two galaxies and was covered with 
complete sampling. North of M31, we made a series of pointed 
observations at locations marked with boxed crosses. The GBT 
angular resolution is 9', but the boxed crosses are 15 x in size to 
make them easier to see. The boxed circle in the lower right shows 
the resolution of BT04. 


1'6 in right ascension. At the end of the 6° strip the 
telescope was stepped in declination by 3!6 and the scan¬ 
ning direction was reversed. This procedure covered the 
area with Nyquist sampling in declination, and finer sam¬ 
pling i n the moving co ordinate to prevent beam broad¬ 
ening (Mangum et al. 2007). Regions at eastern edge 
of each row - areas identified as having no emission by 
BT04 and confirmed from our observations - were used 
as reference positions. Spectra from the eastern-most 
16' of each strip were averaged and supplied the refer¬ 
ence spectrum for the rest of the spectra in the strip. As 
the reference spectrum had an integration time of 30 sec¬ 
onds, this greatly reduced the noise in the difference spec¬ 
trum. The field was observed over and over again until 
the desired noise level was reached. In all, the M31-M33 
field was observed for about 400 hours with an average 
time per GBT beam of 46 min. Our observing procedure 
could cancel some emission, but only if it had a nearly 
constant amplitude and Vlsr over scales of many de¬ 
grees. We see no evidence of this in the data. Th e re are 
three positions in our map where Lockman et al. (2012) 
made frequency-switched detections of HI emission as¬ 
sociated with the clouds. Those data are in reasonable 
agreement with our current measurements. At another 
position where only a small upper limit was reported, we 
likewise see no emission. This indicates that our position- 
switched technique has not cancelled significant amounts 
of HI . In addition, the current data are consistent with 
the sensitive GBT spectrum reported in BT04, but note 
that the declination reported in BT04 has a typograph- 
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ical error: the correct position of the GBT spectrum is 
J2000 = 01 h 20 m 29" + 37°22 / 33 // . 

Spectra were corrected for atmospheric attenuation 
and a second-order polynomial was fit to emission-free 
channels to provide statistics for a quick check on data 
quality. In general, instrumental baselines were exellent 
and modeled well by a 2nd or 3rd order polynomial. A 
small fraction of the spectra (4%) was rejected for hav¬ 
ing poor instrumental baselines caused mainly by radio 
frequency interference or temporary instrumental effects. 
The spectra were smothed to an effective velocity resolu¬ 
tion of 5.15 km s -1 and gridded into a cube using AIPS * 2 
with a pixel spacing of 1.'75 usin g a sp herical Besse l in¬ 
terpolation function following Mangum et al. (2007). A 
third-order baseline was removed from each spectrum in 
the cube. The procedure of subtracting a nearby ref- 
erence position should effectively remove stray radiation 
(Boothroyd et al. 2011) hence no further corrections were 
applied. The noise in the final cube varies slightly with 
position, with a typical value, in brightness temperature, 
of <7 t = 3.45 mK in a 5.15 km s -1 channel. This gives a 
5cr limit on Nhi of 3.9 x 10 17 for a 30 km s -1 FWHM 
line. The equivalent H I mass limit is rsj 10 4 M 0 , as¬ 
suming a distance of 800 kpc. We note that previous 


surveys aroun d galaxy groups t ypically have m ass lim- 
its > 10 5 M^ (|Auld et al.||2006|) or 10 6 (e.g. |Zwaan 

2001 1 |Pisano et al.||2007| |Chynoweth et al.)|2008| . 


2.2. MSI North Observations 

Observations were also made at 18 positions north of 
M31 to investigate the nature of the very faint 21cm 
emission detected by BT04 in this area. The observed 
positions are shown in Figure [l] Here, because reference 
positions were not readily available, the data were taken 
by frequency switching between the 21cm rest frequency 
and a band 4 MHz (844 km s -1 ) away within the 12.5 
MHz band of the GBT Spectrometer. Data were cali¬ 
brated and corrected fo r stray radiation as described by 
Boothroyd et al. (2011), and a third or fourth order poly¬ 
nomial was fit to emission-free velocities. The data were 
then smoothed from 0.3 km s -1 velocity resolution to 
1.3 km s -1 . The typical on-source integration time for 
each pointing is 76 minutes and the median noise in a 1.3 
km s -1 channel is 3.7 mK. This gives a 5 a limit on Nhi 
of 1.6 x 10 17 cm 2 for a spectral line with a full width 
at half maximum (FWHM) of 30 km s -1 . 


3. NEUTRAL HYDROGEN BETWEEN M31 AND M33 

An integrated intensity map of the spectra summed 
over —359 < Vlsr < —187 km s -1 and converted to 

column density, Nhi, is presented in Figure [2] This 
range encompasses all detected HI emission not associ¬ 

ated with the Milky Way. Emission from the disks of 
M31 and M33 overlaps that from the Milky Way at some 
Vlsr > —150 km s -1 , but these velocities are not in 
the range studied here. Nhi is calculated under the as¬ 
sumption that the emission is optically thin, an excellent 
assumption for lines with Tb < 100 mK. The rms noise 
in column density is 1.9 x 10 17 cm -2 (lcr); contours are 
drawn in multiples of 5 x 10 17 cm -2 . The HI emission in 
this field is dominated by discrete clouds, some of which 


2 Developed by the National Radio Astronomy Observatory: 
http:/ / www.aips.nrao.edu/index.shtml. 


are resolved by the GBT 9H beam. Six of these clouds 
were detected at full angular resolution in Paper I, the 
other three appear in the more sensitive data presented 
here. Spectra toward the peak Nhi of each cloud are 
shown in Figure [3j 

A Gaussian function was fit to the 21cm spectrum of 
each cloud at the location of its peak Nhi , and the com¬ 
ponents, along with the location of the peak, are given in 
Table [U Errors are la from the Gaussian fit. Clouds are 
numbered for convenience; the identification is usually 
not the same as in Paper I. 

Table [2] gives derived properties of the clouds. The ve¬ 
locity with respect to the Local Group (Vlgsr) was de¬ 
termined using the calculator given in NED 3 . The quan¬ 
tity Diam is the maximum cloud extent measured down 
to the 3cr noise level while (Diam} is the diameter of a cir¬ 
cle with an area equal to that of the cloud. The difference 
between these quantities is a measure of the elongation 
of a cloud. The “average FWHM” of the 21cm line aver¬ 
aged over the entire cloud is given in Col. 5. Cloud 8 was 
not completely mapped so some of its quantities are lim¬ 
its. There is an additional emission feature at the very 
north-west edge of the field, J2000 = 01 /l 00 m 24 s +37°53' 
with Tl ~ 10 mK at Vlsr = —264 km s -1 and a peak 
Nhi< 10 18 cm -2 , but we cannot characterize it further. 

The HI mass in Col. 7 was calculated by integrating 
over velocities relevant to each cloud, then summing over 
an area around the cloud. The mass assumes that the 
emission is optically thin. Errors on the mass are derived 
from the noise in emission-free channels over the area of 
each cloud. 

A two-dimensional Gaussian was fit to each cloud 
yielding a major and minor axis size. For this, the two 
parts of Cloud 3 were treated separately. Major and 
minor axis radii from the Gaussian fitting were then de¬ 
convolved to produce an estimate of the true angular size 
given the 9'1 beam of the GBT, which is approximately 

Gaussian: r true = (r 2 bs — 4.'55 2 ) 7 . The minor axis of 
Clouds 1, 3b and 4 were unresolved within the errors, 
and for these we assume an intrinsic radius of 1', equiv¬ 
alent to 233 pc at the assumed 800 kpc distance of the 
clouds. Quantities derived from this adopted radius are 
consequently uncertain, and are marked with a colon (:) 
in Table [2| The square root of the product of the decon¬ 
volved major and minor axis radii is given in Col. 6 of 
Table Inland called ri/ 2 - This is the average radius within 
whichnalf the HI mass is contained as estimated from 
the deconvolved Gaussian fit. 

Using this radius we calculate a dynamical (or virial) 
mass - the total mass (from whatever source) needed 
to bind a cloud of radius y 1 / 2 that has a given velocity 
dispersion. We note that the definition of the dynamical 
or virial or total mass derived from the size and velocity 
structure of an object differs considerably from author to 
author, and an exact determination requires information 
that we do not possess, such as the density structure 
within a cloud. For simplicity, and to allow comparisons 
with other measurements (see ^6| we adopt the following 

3 The NASA/IPAC Extragalactic Database is operated by the 
Jet Propulsion Laboratory, California Institute of Technology, un¬ 
der contract with the National Aeronautics and Space Administra¬ 
tion. 
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Figure 2. Integrated HI column density map of the M31-M33 clouds over —359 < Vlsr < —187 km s -1 . The contours are at 
4, 6 and 10 times increments of 5xl0 17 cm -2 . The circle in the lower left shows the angular resolution of the GBT. 
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Table 1 

Clouds Detected Between M31 and M33 


Cloud 

(1) 

J2000 

(hh:mm:ss dd:mm:ss) 

(2) 

Tl 

(mK) 

(3) 

FWHM 
(km s -1 ) 

(4) 

Vlsr 

(km s -1 ) 

(5) 

Nhi 

(10 18 cm- 2 ) 
(6) 

1 

01:24:41.6 +37:24:00 

44.0+2.4 

25.5+1.6 

-297.8 + 0.7 

2.2 ±0.1 

2 

01:23:21.7 +37:18:45 

10.9+2.2 

39.3+9.3 

-223.0 + 4.0 

0.8 ±0.1 

3 

01:20:51.8 +37:15:15 

63.2+1.8 

27.2+0.9 

-237.3 + 0.4 

3.3 ±0.1 

4 

01:19:15.8 +37:29:15 

28.5+2.1 

38.1+3.3 

-228.2 + 1.4 

2.2 ±0.1 

5 

01:16:53.7 +36:49:00 

18.3+1.7 

26.2+3.0 

-308.7+ 1.2 

0.9 ±0.1 

6 

01:08:29.6 +37:45:00 

81.3+3.6 

32.0+1.9 

-278.6 + 0.6 

5.0 ±0.1 

7 

01:05:00.3 +36:21:00 

34.7+3.0 

31.4+3.5 

-210.2 + 1.2 

2.1 ±0.2 

8 

01:03:01.8 +36:00:00 

35.0+2.6 

36.8+3.4 

-281.4 + 1.3 

2.5 ±0.2 

9 

01:01:24.6 +36:12:15 

25.4+3.2 

19.2+2.9 

-341.0 + 1.2 

0.9 ±0.1 


Table 2 

Derived Cloud Properties a 

Cloud 

(1) 

Vlgsr 
(km s -1 ) 

(2) 

Diam b 

(kpc) 

(3) 

(Diam) c 

(kpc) 

(4) 

(FWHM) d 
(km s -1 ) 

(5) 

D/2 e 

(kpc) 

(6) 

Mhi 

(10 4 M 0 ) 

(7) 

Mrlyn 

(io 8 Me) 

(8) 

P s 

(kpc) 

(9) 

1 

-61 

4.4 

3.2 

22.2 + 0.6 

0.38: 

12.7 + 0.2 

0.5: 

126 

2 

+14 

4.4 

2.8 

28.0 + 3.5 

0.75 

4.5 ±0.2 

1.5 

123 

3 

+1 

7.2 

5.3 

27.6 + 0.8 


33.0 + 0.3 


118 

3a 




28.0 + 1.0 

0.82 

22.7 + 0.2 

1.6 


3b 




27.1 ± 1.3 

0.34: 

10.2 ±0.2 

0.6: 


4 

+12 

3.7 

3.1 

26.5 + 1.5 

0.41: 

8.6 ±0.2 

0.7: 

112 

5 

-69 

7.0 

3.8 

21.8 + 2.5 

1.13 

7.8 ±0.2 

1.3 

109 

6 

-32 

7.2 

4.6 

33.6 + 1.3 

0.78 

39.2 ±0.3 

2.2 

87 

7 

+36 

4.9 

3.9 

27.4 ±2.1 

0.71 

12.6 ±0.4 

1.3 

92 

8 

-35 

> 3.5 

> 3.2 

29.6 + 2.1 


> 11.6 


91 

9 

-94 

3.5 

3.2 

20.5 + 2.3 

0.96 

8.7 ±0.2 

1.0 

88 


a For an assumed distance of 800 kpc. 
b Maximum cloud extent. 

c Diameter of a circle with an area equal to that of the cloud. 
d Column density weighted average FWHM over the entire cloud. 
e Square root of the product of the major and minor axis radii. 
f From Equation]^ 
g Projected distance from M31. 
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Figure 3. Spectra taken at the peak Nhi of the nine clouds de¬ 
tected in the M31-M33 field. The vertical axis is 21cm brightness 
temperature and the horizontal axis is Vlsr- 


(Binney & Tremaine ^008): 


Mdyn ^ 


( 1 ) 


where G is the gravitational constant and cr v is the 3- 
dimensional velocity dispersion of the cloud. As cr v is 
not a measurable quantity, we use the FWHM from the 
measured average spectrum, and assume isotropy to cal¬ 
culate 


Mdy 




-yn \ _ 


m 2.5 x 10^ 


/r 1/2 \ / FWHM V 
ykpc J \ km s _1 J 


( 2 ) 


where ii/ 2 is in kpc, FWHM in km s -1 and Md yn is in 
solar masses. Dynamical masses are given in Col. 8. In 
all cases they exceed the observed HI mass by a factor of 
~ 10 3 . Col. 9 gives p, the projected distance from M31 
assuming that the clouds are at a distance of 800 kpc. 
The clouds lie > 90 kpc away from M31’s center, farther 
than the distances of the known HVCs around M 31, all 
of which have p < 50 kpc ( Westmeier et al.]|2008 ). 


3.1. The Velocity Range 

In the field of Figure [2j we find no HI emission at 
velocities more negative than that associated with cloud 
9, i.e., nothing at Vlsr < —370 km s -1 and no emission 
at the more positive velocities between —190 < Vlsr 


— 160 km s 


-l 


At still more positive Vlsr, there is a 


band of HI extending from the south center of the field 
to the north-west that is quite bright by our standards 
(T b > 0.6 K) at -150 < V LS r < -100 km s” 1 . This 
appears to connect smoothly to Milky Way emission at 
more positive velocity and we will not consider it further 
here. 


Lehne r et al.| (|2015) argue that the CGM of M31 should 
be defined as having —300 < Vlsr ^ —150 km s -1 , 
that velocities more negative than this may arise in the 
Magellanic Stream and velocities more positive in the 
halo of the Milky Way. Our Cloud 9 has a peak Nhi at 
Vlsr = —341 km s _1 , but is in no way unusual in its 
size, mass, linewidth, or projected distance from M31, 
so we assume th at it is part of the cloud population. 
Westmeier et al. (2008) have detected 21cm emission that 


they attribute to M31 HVCs only ^ 2° to the northwest 
of the Figure [2] field at velocities as negative as Vlsr ~ 
—500 km s _1 . We would have easily detected similar 
emission in the Figure [2] field. Apparently the M31 HVC 
population does not extend over this area. This will be 
discussed further in © 

3.2. The Total Neutral Hydrogen Mass 

Figure [2] shows some evidence for HI emission outside 
of the clouds listed in Table [lj but this is relatively small 
and ususally concentrated in discrete regions. Assuming 
that all the HI we detect is at a distance of 800 kpc, 
the sum of the HI mass of all nine clouds is 1.4x 10 6 
Mq, while integration over the entire field yields 1.6 x 10 6 
M©. Thus virtually all of the neutral gas we measure is 
contained in the nine clouds with about half in just the 
two largest clouds, numbers 3 and 6. 

To compare our data with those of BT04, we match the 
angular resolution of that survey by convolving the inte¬ 
grated spectra shown in Figure [2] to an ang ular resolu¬ 
tion of 49' using the WSRT beamlrom Popping_& Braun 

:e [4] While the maps 


(2008). Results are shown in Figure 


have many similarities there are significant differences. 
The most striking is near J2000 = 01 /l 13 m 30 s ,+37 0 24', 
where the BT04 data show emission that does not ap¬ 
pear in the GBT data. The BT04 survey also contains 
emission at the south-central part of our field, some 
of which may arise from a cloud just off our map at 
01^10 m ,+35°30 / convolved with the much larger BT04 
beam. In some areas we find good agreement between 
the two surveys - the northwest and northeast corners, 
for example - but overall BT04 reports an HI mass 
2.5 x 10 6 M© for the area, while we find only 63% of 
this amount. 

We believe that the cause of this discrepancy lies 
mostly, though not entirely, in the choice of velocity 
range of integration for the BT04 survey. No emission is 
detected with the GBT at velocities more negative than 
—370 km s _1 . If the GBT data are integrated over all ve¬ 
locities < —150 km s -1 , the total mass remains 1.6 x 10 6 
M©, but if the upper limit is taken to be —140 km s -1 , 
only 10 km s -1 more positive, the total HI mass doubles. 
We discuss above the reasons for not including emission 
at Vlsr > —150 km s -1 in the census of L ocal Gro up 
gas, a s it likely arises in the Milky Way (Lehner et al. 


(2015) conclude this as well). It is plausible that some of 
this emission is present in the BT04 map, especially as 
those spectra had a rather coarse velocity resolution of 17 
km s -1 . This explanation does not account for the dis¬ 
crepancy with the BT04 feature at 01 /l 13 m 30 s + 37 0 24'. 
The GBT spectrum in this direction is entirely consistent 
with noise at Vlsr < —145 km s -1 . 

In Paper I we reported that there was « 1 x 10 6 M© 
of HI in the GBT data that was not associated with 
discrete clouds. We now believe that this conclusion is 
incorrect and resulted from a very small systematic base¬ 
line error with an amplitude of only a few mK in the 
preliminary GBT data. In the current data there can 
be no more than 0.2 x 10 6 M© of HI outside the nine 
clouds. We do not believe that our data reduction proce¬ 
dure is artifically suppressing real emission; instrumental 
baseline fitting removes only low-order polynomials pre¬ 
serving lines of normal velocity width, and the position¬ 
switching observing technique preserves structure on an- 
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Figure 4. Top panel: integral column density of the M31-M33 
field at full GBT angular resolution; contours are every 5 X 10 17 
cm -2 begining at 5 X 10 17 cm -2 and the color scale runs from 
zero to 4 x 10 1 ® cm -2 . Middle panel: GBT data smoothed to the 
49' angula r resolution of the BT04 m easurements using the WSRT 
beam from Pop ping fc Braun] (2008 ). For this and the lower panel, 
the color scale runs from zero to 4 x 10 17 cm -2 , and the contours 
are every 10 17 cm -2 beginning at 10 17 cm -2 . Lower panel: BT04 
integral column density map at 49' resolution. Regions without 
detectable HI are blanked. While many of the general features of 
our GBT data agree with those from BT04, the BT04 map contains 
regions of emission not found in the GBT data. 

gular scales <5°. If there were emission in the reference 
regions beyond the edge of the map, we would see it as 
negative features in the map. As no such features exist, 
we have confidence in our estimates of Mhi . 



01 h 28 m 24 m 20 m 16 m 12 m 08 m 04 m 00 m 


J2000 Right Ascension 

Figure 5. Average Vlsr over the field, calculated for channels 
between —359 < Vlsr < —187 km s -1 that have > 4cr = 0.015 
K. Although velocities can vary significantly from cloud to cloud, 
the mean velocity within a cloud is relatively constant. The cir¬ 
cle in the lower left shows the angular resolution of the GBT. It 
is apparent that the clouds are separated not only spatially but 
kinematically. 


Cloud 3 

i -i-1-i-i-i-1-1-1-r 



01 h 22 m 40 8 00 s 21 m 20* 20 m 40 8 00* 


J2000 Right Ascension 

Figure 6. Average Vlsr pixel-by-pixel for Cloud 3 together with 
contours of Nhi • Averages were calculated only for channels with 
Ts > 4cr = 15 mK. Contours are at 1, 2, 3, 4, and 5 times in¬ 
crements of 5 X 10 17 cm -2 . Unrelated emission from Cloud 4 
appears in the upper right. The circle in the lower left shows the 
angular resolution of the GBT. Both cloud components show a 10 
km s -1 velocity gradient between their center and edge. 


3.3. Velocity and Linewidth 

The average velocity of the HI emission is shown pixel- 
by-pixel in Figure [5] Contours are the same as in Fig¬ 
ure [2] Here, to give an indication of the overall velocity 
pattern of the clouds the velocity-weighted values of the 
brightness temperature were calculated over the entire 
range of detectable emission: —359 to —187 km s -1 , 
clipping the data at the approximate 4cr noise level of 
15 mK. The clouds span a range of velocity between 
—341 < Vlsr < —210 km s -1 . In our data, there are 
only two clouds that have significant internal velocity 
structure. These are shown in Figures [6] and [7] 

Cloud 3 consists of two components with nearly iden¬ 
tical velocities at their peak, but overall there is a 10 
km s -1 gradient from —240 km s -1 between the two 
cloud components to —233 km s -1 at the northern edge 
and —230 km s -1 at the southern edge. Cloud 6 has 
a head-tail appearance with, once again, velocities most 
negative at the center and rising by « 10 km s -1 at the 
northern and southern edges. 

The mean linewidth at the peak Nhi of the nine clouds 
is 27.9 d= 4.3 km s -1 , essentially identical to the me¬ 
dian. The lines are well-fit by single Gaussians. Unlike 
Galactic HVCs and CHVCS, there is no evidence for two 



01 h 10 m 00 8 09 m 20 8 08 m 40 8 00 8 07 m 40 8 

J2000 Right Ascension 


Figure 7. Average Vlsr pixel-by-pixel for Cloud 6 together with 
contours of Nhi- Averages were calculated only for channels 
with Tl > 4cr = 15 mK. Contours are at 1, 2, 3, 4, 5, 7.5 and 
10 times increments of 5 X 10 17 cm -2 . The circle in the upper left 
shows the angular resolution of the GBT. This cloud shows a 10 
km s -1 velocity gradient from the center to its edges. 
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Figure 8. Expanded view of the BT04 survey HI data from Fig- 
ure^to show the locations of the GBT measurements to the north¬ 
west of M31. The GBT angular resolution is 9', but the symbols 
are 15' across for easier identification. Green boxes mark positions 
with detected HI emission at Vlsr < —150 km s -1 . Red circles 
show positions without a detection. The circle in the lower right 
shows the BT04 angular resolution. 


components in any of the lines as would be e xpected if 


they have a two-phase temperature structure (Sternberg 
et al. 2002| ). The FWHM does not vary much across 


an individual cloud. The FWHM averaged over the en¬ 
tire cloud (Table [2k and toward the peak N H i (Table [IJ 
are quite similar. The range of linewidths is 19.2 - 321) 
km s -1 , much smaller than the range observed in 17 M31 
HVCs (at e ssentially identical ang ular resolution) of 11 - 
71 km s -1 ( |Westmeier et al.||2008 ). 

Figure [5] re-enforces the impression from Figure [2] that 
the individual clouds are independent entities, as nearby 
clouds can have quite different velocities. Clouds 7 and 
9, for example, differ by 120 km s -1 yet are separated by 
< 1° on the sky (^ 10 kpc at 800 kpc distance). 


4. NEUTRAL HYDROGEN TO THE NORTH OF M31 

Results from the measurements to the north of M31 are 
given in Table [3j and illustrated in Figure [8] as symbols 
on a map of the BT04 survey. We include only emission 
with Vlsr < —150 km s -1 . Line properties were de¬ 
rived from a Gaussian fit, and errors are la. Values of 
Nhi come also from the Gaussian fit; here errors reflect 
both the noise and an assumed equal contribution from 
baseline uncertainties. The median 5 a sensitivity to a 30 
km s -1 FWHM line isl.6xl0 17 cm -2 . Because we have 
only incomplete sampling, it is not possible to delineate 
objects and calculate a mass or size. The values of Nhi 
should be understood as random samples of the medium 
and not peak values, in contrast to the values given in 
Table □ for the M31-M33 field. 

Our results are rather puzzling in view of the BT04 



-450 -400 -350 -300 -250 -200 -150 -100 -50 

V|_sr (k m s 5 

Figure 9. [a] Spectrum from a deep GBT pointing toward J2000 
= 00^15 m + 46°00 / . The component at —238 km s -1 has a velocity 
similar to M31’s HVCs while the component at —388 km s —1 
is closer to the systemic velocity of M31. [b] Spectrum from a 

pointing toward J2000 = 00^20 m + 47° 00 7 . The spectral line lies 
~ 50 km s —1 from Milky Way emission. 

data. We detect no emission at S = +45°, at only 
one position at S = +46°, and yet have four detec¬ 
tions at S = +47°, two of which lie outside the lowest 
BT04 contour. The first conclusion must be that the 
HI is much patchier than would be inferred from BT04, 
which does not give a good representation of the HI at 
Vlsr < —150 km s _1 . In this sense the M31 north mea¬ 
surements complement the conclusion from the M31-M33 
field (Figure El). In the GBT data, bright lines iden¬ 
tified with Milky Way emission are regularly found at 
Vlsr ~ —100 km s -1 . It is possible that the BT04 map 
includes some of this material, which might account for 
the discrepancy with the GBT measurements. Our de¬ 
tection of HI outside the BT04 contours at 5 = +47° 
would imply that we are seeing very small angular sized 
features that suffer from beam dilution in the BT04 mea¬ 
surements. This was also a conclusion from Paper I. 

The position at 00H5 m , +46 o 00' has two 21cm line 
components separated by 150 km s -1 . Figure [9^i shows 
this spectrum as well as the spectrum at 00 /l 20 m , +47°00' 
®>). Both illustrate the emission near —100 km s -1 that 
we attribute to a component of the Milky Way. 


5. KINEMATICS OF THE EMISSION 

Velocities of all GBT detections are shown in Figure [To| 
as the Local Group Standard of Rest velocity (Vlgsr) 


versus angular distance from M31. We 
clouds of Figure [2] and the detections 
M31. For comparison, the known HV 
and M33 are indicated with red circles ( 

include both the 
to the north of 
r Cs around M31 
Westmeier et al. 

2008 Grossi et al. 

2008; 

Putman et a 

2009) and the 

systemic velocities of 1VI31 and 1VI33 witJ 

ti blue rectangles 


(Karachentsev & Makarov [1996). As discussed in Pa¬ 
per I, the clouds between M31 and M33 lie at velocities 
similar to the systemic velocity of both galaxies, and at 
a larger distance from either galaxy than their HVCs. 
This is true also of the two newly-identified clouds be¬ 
tween M31 and M33. Most of the emission to the north 
of M31 appears consistent with arising in an extension of 
the M31 H VC p o pulation . It is in the velocity range as¬ 
sociated by Lehner et al. (2015) with the circumgalactic 
medium of MUT 





























































Table 3 

Observations North of M31 


J2000 

Galactic 

O t 

T l 

Vlsr 

FWHM 

Nhi 

Vlgsr 

P 

(hh:mm dd:mm) 

1°, b° 

(mK) 

(mK) 

(km s -1 ) 

(km s -1 ) 

(10 17 cm 2 ) 

(km s -1 ) 

(kpc) 

(i) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

00:10 +46:00 


3.8 







00:10 +47:00 

115.6, -15.3 

4.4 

27.5 + 1.3 

-176 + 1 

25.9 ± 1.4 

13.8 + 1.0 

106 

115 

00:15 +46:00 

116.3, -16.4 

3.5 

15.6 + 1.2 

-238 ± 1 

16.0 ± 1.4 

4.8 ±0.6 

43 

96 



3.5 

22.4 + 1.2 

-388 ± 1 

23.6 + 1.4 

10.2 + 0.8 

-109 


00:15 +47:00 


3.3 







00:20 +45:00 


6.4 







00:20 +47:00 

117.4, -15.5 

3.3 

30.4 + 1.0 

-166 + 1 

26.1 ± 1.0 

15.4 + 0.7 

112 

98 

00:20 +48:00 


3.6 







00:25 +45:00 


3.8 







00:25 +47:00 


3.7 







00:30 +45:00 


4.1 







00:30 +46:00 


4.4 







00:30 +47:00 


4.5 







00:35 +45:00 


4.1 







00:35 +47:00 


3.7 







00:40 +45:00 


3.3 







00:40 +46:00 


3.8 







00:40 +47:00 

120.9, -15.8 

2.9 

18.3 + 0.7 

-182 + 1 

25.1 ± 1.1 

8.9 ±0.6 

85 

81 

00:45 +47:00 

121.8, -15.9 

3.7 

18.8 + 0.8 

-183 + 1 

42.5+ 2.1 

15.5 ±1.0 

86 

80 



Figure 10. Velocity with respect to the Local Group Standard of 
Rest (Vlgsr) versus angular distance from M31, where directions 
toward M33 are taken to be positive. The blue squares are M31 and 
M33, with the red dots being the high velocity cloud populations 
of each galaxy. The black crosses are the clouds from the M31- 
M33 map with vertical lines marking the map limits. The green 
diamonds are detections from the GBT pointings north of M31. 
The pink cr oss marks the UV absorption line measurements of 
|Koch et al.| ( |2015| > towards a backgrou nd AGN. The MS curve i s 
the estimated position of Stream SO from |Nidever et ak] (|2010), with 
=b 30 km s -1 spread to account for the velocity width of the Stream 
and our approximations. In most cases the Stream lies hundreds 
of km s -1 away from the emission we detect. The rotation of M31 
extends from -353 to +248 km s~* in Yrngp. and for M33 from 
-105 to +158 km s -1 ( |Braun et al.|2009| Putman et al.|2009| . 


Koch et al. (2015) report detection of gas possibly as¬ 
sociated the the 1V131 Stellar Stream through UV absorp¬ 
tion line measurements against background AGN. Their 
most secure detection is in lines of Si, C, and O against 
Q0045+3926 at a location about 2° from M31 in the di¬ 
rection of M33. They find one system at Vlsr = —370 
km s _1 , or Vlgsr = — 111 km s -1 , which is marked with 
a pink cross in Figure [To| It lies in the region occupied 
by M31 HVCs. 

The second HI component toward OOMS 7 ™, +46°00' 
(Figure [ 9 J 1 ) at V L sr = -388 km s -1 (Vlgsr = -109 


km s -1 ) has a velocity that might be ass ociated with 
the Magellanic Stream (Lehner et al. 2015) though this 
direction lies nearly 2 0° from the ax i s of th e Magellanic 
Stream as defined by Nidever et al. (2010), and not far 
from a direction t hat does not sh ow U V absorption aris¬ 
ing in the Stream ( Fox et al. 2014). As our measurements 
are only over a sparse grid in tins area, we do not know 
the size or mass of any cloud associated with this fea¬ 
ture. One of the M31 HVCs, the Davies Cloud, has also 
been suspected of being part of the Stream because its 
mass would be an order of magnitude larger than the 


other M31 HVCs were it at the distance of M31 (Davies 
1975; Westmeier et al.|2 008). The Davies Cloud does not 


have extraordinary kinematics and does not occupy an 
unusual location in Figure 
its HI mass is similar to t 


10 At the distance of M31 


hat of the Milky Way HVC 


Complex C (Wakker et al. 2007a). We see no reason to 
exclude it from the 1V131 HVC po pulation. 

Because Nidever et al. (2010) have shown that the 
Magellanic Stream (MS) extends into the BT04 area, we 
have tried to estimate where the MS would lie in Fig¬ 
ure m Five HI Streams wer e identified in the BT04 

Figure 4), labeled S0-S4. 


data from Nidever et al. 


(2010 


We take Stream SO, the one""closest to our pointings, to 
lie along a line of constant RA at 23 /l 45 m and choose a 
declination range from +20 to +50 degrees. 

To approximate the LGSR velocities of Steam SO along 
this line, we first transform the positi ons to Magellanic 
Longitude (Lms> [Nidever et al. 2008) and use the ex¬ 
tended Vlsr vs. Lm s fidncTal curve from Figure 7b of 


Nidever et al. ([2010), converting to Vlgsr- The result¬ 
ing angle-velocity for Stream SO is shown as a dashed 
curve in Figure [l0| To account for the spread in velocity 
of the Stream, and for our approximations, we over-plot 
a shaded region with bounds of ± 30 km s -1 from the 
curve. 

The point of closest approach of SO to M31 is at an 
angular distance of 10.3° with Vlgsr = -187 km s -1 . 
Most of our detections to the northwest of M31 lie at least 
225 km s -1 from the average velocity of SO (Vlgsr ~ -175 
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km s- 1 ). The line at 00M5 m + 46°00', Vlgsr = -109 
km s _1 is the only one within 5° and 50 km s _1 of SO, 
so excepting this, it seems unlikely that the emission we 
detect arises in the Magellanic Stream. 

6. DISCUSSION: THE ORIGIN OF THE CLOUDS BETWEEN 
M31 AND M33 

In the 12 square-degree area between M31 and M33 ob¬ 
served with the GBT, the faint HI discovered by BT04 
is resolved into at least nine discrete clouds. Here we 
consider possibilities for the origin of the clouds. 


6.1. Products of Galaxy Interactions 

Although BT04 surveyed a large area, they found HI 
at relevant velocities only near a line joining M31 and 
M33 on the sky (Figure fTI) , and not widespread in the Lo¬ 
cal Group. This naturally leads to models where the HI 
tracks past tidal interactions, as in t he M81 group (Yun 
et al. || 19941 |Chynoweth et al. j2008|) . |Bekki| (|2008 ) mod¬ 
eled the BT04 results as arising from a close encounter 
between M31 and M33 some 4-8 Gyr ago, and others 
have also considered an interaction a few Gyr ago as an 
explanation for features in M33’s extended HI distribu¬ 
tion and the stellar structures in M31 and M33 d Putman 


et al.||2009| |McConnachie et al.||2009| |Lewis et al.||2013p 

The Magellanic Stream is a good, nearby template of ai 


good, near by template of an 


ongoing galaxy intera ction (e.g., Besla et al. 2010 Sta- 
nimirovic et al. 2010) but its age is quite uncerta in, m 
the range of 0.3 - 2.0 Gyr ( jD’Onghia fe Fox |2015|). The 


HI mass and average Nhi is gener ally higher than we 
observe between M3 1 and M 33 (Stanimirovic et al.[2008; 
Nidever et al.||2008[ |2010|). It is worth noting that as 
gaseous tidal debris ages it should disperse causing its 
Nhi to decrease, so the M31-M33 clouds may represent 
material that is substantially older than is seen in the 
Magellanic Stream. 

A recent study of the evolution of the Local Group has 
now cast doubt on the possiblity that we are observing 
the relic of a past interaction between M31 and M33, as 
it indicates that in the last ~ 12 Gyr, M33 has not been 
closer to M31 than it is now (|Shaya & Tully||2013|). The 
same calculations though, suggest that the dwarf galaxies 
And II and And XV could have interacted with M31 over 
the last 0.7 Gyr creating an HI extension toward M33. 
Since Local Group dwarfs with detectable HI typically 
have Mni/Ly ~1 (McConnachie| 2012 Spekkens et al. 
2014|), we can estimate how much H1 And 11 and And 
XV could have contributed to the M31 CGM. Based on 


the My from McConnachie (2012), their total HI mass 
would be ^ 10 Y M©, which is comparable to the total HI 
mass in the entire filamentary structure between M31 
and M33. It seems unlikely, however, that the encounter 
would leave the nine clouds we detect, each apparently 
coherent, but spread over a projected distance of 70 kpc 
with a spread in Vlgsr of 130 km s _1 . It is clear that 
more detailed modelling of the Local Group, including 
its gas, is needed to address these possibilities. 


6.2. Dwarf Galaxies 

The larger M31-M33 clouds have an HI mass similar to 
that of some dwarf galaxies, but the clouds are not likely 
to be associated with stellar systems. Table [4] compares 
properties of the most massive cloud with dwarf galaxies 


of a similar HI mass. The dwarf spheroidal And XII is 
included as it is one of the faintest of the known M31 
satellites, but was d etected easily in a recent search for 
even fainter systems (Martin et al |2013). To allow for ac¬ 
curate comparisons we calculate"dynamical masses from 
the measured quantities in Table [4j using Eq. [2] 

Local Group dwarfs with detectable HI ty pically have 


a stel lar mass similar to their gas mass (Spekkens et al. 


2014). For Cloud 6 this would imply My < —8, and its 


stars would certainly have been detected already. |Martin 


et al. (2013) report a possible stellar feature at the Cloud 


6 positionTout is considerably fainter than And XII, and 
if it is actually associated with the Cloud, would imply 
Mri/M* > 10. The dynamical mass of Cloud 6 is also 
about an order of magnitude higher than that of Leo P, 
even though Leo P has ^ 10 times more baryonic mass. 
Using the relationship bet ween total mas s and My from 
the Local Group data of McConnachie (2012), Cloud 
6 - indeed most of the clouds detected in HI between 
M31 and M33 - should have a stellar counterpart with 
My < —13, whereas the surveys of stars near M31 sug¬ 
gest My > — 6. If Cloud 6 is a galaxy it is one with 
rather extreme properties. It will be interesting to see if 
the possible a ssociat ion wit h a sli ght stellar overdensity 
discovered by Martin et al. (2013) reveals a real stellar 
component of this cloud. The dwarf galaxy Leo P has 
a small rotational velocit y, 15 km s -1 , comp arable 
to its velocity dispersion ( Bernstein-Cooper et al.||2014 ). 
As shown in Figure [7], Cloud 6 shows a slight gradient in 
Vlsr from center to edge; at our angular resolution any 
rotational component is comparable to or smaller than 
its velocity dispersion. 

Another reason why the M31-M33 clouds are not likely 
to be galaxies is that the stellar satellites of M31 and the 
Milky Way that lie closer than ss 300 kpc to the parent 
galaxy are extremely deficient in HI , with mass limits 
typic ally well below 10 4 M©. In some cases Mht ^ 100 


M© (Grcevich & Putman 2009 Spekkens et al. 2014 
R. Beaton private communication). The two 1V131 satel- 
lite galaxies that lie nearest the M31-M33 field, And 
II with Vlsr = —187 km s -1 , and And XV with 
Vlsr = —322 km s" 1 , are amo ng those that lack de¬ 
tectable HI emission (Lockman et al. 2012 ] ). Apparently 
M31 and the Milky Way are very efficient at stripp ing gas 
from small s atellites passing through their CGM ([M ayer 


et al. 120061 [Grcevich fe; Putman||2009| |Nickerson et al. 


201 H |Gatto enL|2013). 


Figure |lT[contains the same HI data as Figure [10 but 
instead ofthe M31 and M33 HV Cs, the locat ion and ve¬ 
locity of dwarf galaxies from McConnachie| (2012) are 
shown. There is no obvious connection between the 
dwarf galaxies and the HI clouds as the dwarfs are spread 
over 400 km s -1 while the clouds have an average Vlgsr 
like that of M31 and M33 and a total range of only 130 
km s _1 . 


6.3. High Velocity Clouds and Dark Matter Sub-Halos 

Table [5] gives information that allows us to compare 
the M31-M33 clouds with various populations of HVCs. 
For the M31-M33 clouds, values are the median from Ta¬ 
ble [21 Fo r the M31 HVCs the list ed radius is the median 


found by Westmeier et al. ( 

2005 

) while the FWHM and 

Mh i are the median values from 

Westmeier et al. (2008). 

Values for the Ultra-Compact H 

^Cs (UCHVC) are from 
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Table 4 

Comparison of an M31-M33 Cloud with Dwarf Galaxies 


Object 

(i) 

r l/2 

(kpc) 

(2) 

FWHM b 
(km s -1 ) 

(3) 

M H i 

(M 0 ) 

(4) 

M* 

(M 0 ) 

(5) 

M v 

(mag) 

(6) 

Mrlyn 

(M 0 ) 

(7) 

References 

(8) 

Cloud 6 a 

0.78 

34 

3.9 x 10 5 

— 

— 

2.2 X 10 8 

a 

Leo P 

0.25 

24 

9.5 x 10 5 

5.7 x 10 5 

-9.4 

3.6 x 10 7 

b,c,d 

Leo T 

0.17 

16 

2.8 x 10 5 

1.4 x 10 5 

-8.0 

1.1 x 10 7 

e,f,g,h 

And XII 

0.30 

6 

— 

3.1 x 10 4 

-6.4 

2.8 x 10 6 

i,j 


Referenc es. — (a) This Work; fb)|B ernstein-Cooper et al. (2014}; (c ) [Adams et al. 

7]); (g) |t'aerman et al.| J2U13| >; (h) [<S 


|<AjdT ( p?008| ); (f) |Simon fc Ceha| ( |200 y| >; (g) Faerman et al. ;20l:(p (h) |de Jong et al.| ( f^UU«[ l ; (i) |uollms et al.| 


; (d) [McQuinn et al. 


UsSl»; ( e ) 1 

Arruf ; (j)[M 


Ryan-Weber 

/IcUonnactue 


a For an assumed distance of 800 kpc. 

b Of the HI emission except for And XII where it is from the stars. 



Figure 11. Velocity with respect to the Local Group Standard of 
Rest (Vlgsr) versus angular distance from M31, where directions 
toward M33 are taken to be positive. The blue squares are M31 and 
M33, and the black crosses are the clouds from the M31-M33 map. 
The green diamonds are detections from the GBT pointings north 
of M3 1. The pink trian gles are dwarf galaxies from the compilation 
of|McConnachie| ([2012). There is no apparent connection between 
the dwart galaxies and the M31-M33 HI clouds. 


Adams et al. (2013) scaled to a distance of 800 kpc. 

There is considerable overlap between the physical 
properties of the M31-M33 clouds and the M31 HVCs, 
but much less so between their locations and kinemat¬ 
ics. Many HVCs are located cl ose to th e edge of t he 
M31 disk in intersecting filaments (Westmeier et al. 2005 
. Although median line widths are similar, the M3 
HVCs have widths as low as half and as high as twice 
the FWHM of the M31-M33 clouds; 90% have a FWHM 
in the range 13-67 km s -1 , whereas the entire sample of 
M31-M33 clouds has FWHMs between 19 and 39 km s -1 . 
The HVCs also have a median mass four times larger 
than the clouds. Indeed, there is no cloud in our sample 
that has an HI mass as high as the median HI mass 
of the M31 HVCs. Most importantly, the kinematics of 
the clouds differ significantly from that of the HVCs near 
them (Figure [To]) , and they have a substantially smaller 
spread in velocity (see also Paper I). With only one ex¬ 
ception the detections north of M31 have positions and 
velocities consistent with being part of the M31 HVC 
population while the nine clouds between M31 and M33 
do not. The M31 HVCs in the direction of M33 have 
considerably more negative velocities than the systemic 


velocity of M31. For M31 Vlsgr = —34 km s -1 , while 
its HVCs towards M33 have (Vlgsr) = —166 km s -1 . 
In contrast, the M31-M33 clouds have (Vlgsr) = —25 
km s -1 . The clouds do not appear to be simply an ex¬ 
tension of the M31 HVC population. Still, the physical 
properties displayed in Table [5] suggest that the two pop¬ 
ulations may be formed by similar processes, and that the 
UCHVCs, if at a similar distance, have similar properties 
as well. 


Thilker et al. (2004) and Westmeier et al. (2005) con¬ 
sidered whether the HVCs of 1V131 might be the baryonic 
component of a population of low-mass dark matter ha¬ 
los, an idea with a long history in the study of HV Cs 
(|Oort||1966| |Blitz et al.||1999| iBraun & Burtonll 19991 IdF 


Heij et al. 2002J Mchols et al. 2014). Of particular im- 

portance is that some HVCs (and some dwarf galaxies) 
show a two-component structure in their 21cm lines, sug¬ 
gesting the presence of gas in equilibrium at two temper¬ 
atures. This is a diagnostic of pressure, and was used 
(among other pieces of evidence) to argue that the com¬ 
pact HVCs (CHVCs) must lie at distances < 150 kpc 
from the Milky Way, a nd not at distances of ~ 1 Mpc 
(|Sternberg et al.||2002|). The failure to detect CHVC 
analogs spread throughout other groups of galaxies also 
implies that th ey must be locate d within 90 kpc of indi¬ 
vidual galaxies (Pisano et al. 2007). This result is consis¬ 
tent with HVC detections arouna other galaxies such as 


NGC 891 and N GC 2403 (Oosterloo et al. 2007 Frater^ 
nali et al.|2 002|) and the distance br ackets f or Milky W ay 
HVCs (e.g. jWakker||2001 Wakker et al.||2QQ7b , 2008|). 
The UCHVCs are discussed in depth in I Adams et al. 


HeptlT in |Ad; 

(2013) and|Faerman et al. (20131). The dwarf galaxy Leo 
P was originally classified as a UCHVC on the basis of 


originally classified as a UCHVC on the basis of 
its 21cm characteristics before its stellar component was 
detected and it was discovered to be at a distance of 


1.7 Mpc (Giovanelli et al. 2013 

Rhode et al. 2013 

Me- 

Quinn et al. 2013). The HI properties listed in Table |5 


HVCs and the M31-M33 clouds. There is no evidence 
for a two-phase interstellar medium in the M31 HVCs, 
the M31-M33 clouds, or in the unidentified UCHVCs, 
although it is found in other HVCs and the Magellani c 
Stream ( Kalberla fc Haud|2006 Stanimirovic et al.|2010 ). 


This is important as it implies tn at these ob jects exist in 


regions of low external pressure (Sternberg et akl 2002). 


The observations to date, however, have relatively poor 
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Table 5 

Comparison M31-M33 Clouds with HVCs 


Object 

(i) 

r l/2 

(kpc) 

(2) 

FWHM 
(km s" 1 ) 

(3) 

M H i 

(M 0 ) 

(4) 

M dyn 

(M©) 

(5) 

References 

(6) 

M3f-M33 Clouds 

0.75 

27 

1.2 x 10 5 

1.4 x 10 8 

a 

M31 HVC 

0.52 

26 

4.7 x 10 5 

8.8 x 10 7 

b, c 

UCHVC 

1.16 

23 

1.2 x 10 5 

1.5 x 10® 

d 


References. — (a) This Work; (b) |Westmeier et al.| ( |2005| >; (c) |Westmeier et ah] |2008 ); (d) [Adams et al.| ( |2013| > 


linear resolution and might not be able to detect a cool 
HI phase if present in moderate amounts. 

6.4. Gas in Planes of Satellites or a Dark Matter 
Filament 

Lacking an estimate of the distance to the clouds, it 
is not clear if they are related to the planes of satellite 
galaxies that are now tho ught to be fundamental struc¬ 
tures in the Local Group (Ibata et al.||2013 Conn et al. 
2012 2013). The M31-M3d held studied here lies in the 
part of the sky where M31 satellites are aligned in the 
structure called “Plane 2”. It lies t o the west of M31 and 
extends southwards towards M33 ( Shaya fc Tully||2 013). 
The two M31 satellite galaxies nearest to our held on the 
sky have velociti es within the range o f the HI clouds 
(see Figure 3 of Shaya & Tully ( 2Q13| ) ) and it is an in¬ 
triguing possibility that the M31-M33 clouds are part 
of a larger alignment of matter in this part of the Local 
Group. Although study of these planes is only beginning, 
if they represent large-scale dark matter structures, then 
it is very plausible that they would be accompanied by 
enhancements in gas density. 

BT04 suggested that the HI they discovered origi¬ 
nated from condensation of hot gas in a dark matter 
filament connecting M31 and M33. This is now testable, 
in part, through cosmological simulations that attempt 
to follow the evol utio n of systems like the Loca l Group. 
Nuza et al. (2014) and Scannapieco et al. (2015) have an- 
alyzed the simulation of a group containing two galaxies 
like M31 and the Milky Way, and while M33 is consid¬ 
erably less massive then the Milky Way, the simulation 
might still have some application to the M31-M33 re¬ 
gion. They find that the hot gaseous halos of the M31 
and Milky Way analogs overlap and that they evolve to 
occupy the same filament, leading to an excess of neutral 
gas between them that forms around z ^ 1 and persists 
to the present. Clouds like those observed here might 
then condense in the filament. The current simulations 
do not have the resolution to detect anything as small as 
the M31-M33 clouds, and it is not clear that M31 and 
M33 would have similar overlapping gas halos, but these 
results are encouraging that intra-group gas might be a 
natural feature of systems like the Local Group. 

In view of the evidence that M31 has a massive CGM, 
as discussed in the next section, it will be important to 
determine if the direction of M33 is enhanced in total 
material, as would be suggested by this scenario. 

6.5. Condensations in the MSI Circumgalactic Medium 

While the presence of the clouds between M31 and M33 
is unexpected, the total mass involved is not large com¬ 


pared to the baryonic mass at that location in the halo 
of M31 as determined from recent measurements. The 
observed Nht averaged over our field is only 9 x 10 16 


, —2 


cm " . Lehner et al. (2015) have studied the extended 


CGM of M31 through measurement of UV absorption 
lines against background QSOs out to a projected dis¬ 
tance > 500 kpc, an area that covers not only the M31- 
M33 clouds, but the galaxy M33 as well. They find that 
the total gas mass in the M31 CGM may be > 10 10 
M© with an ionization fraction > 90%. The Lehner et 
al. radial column density profile of Sill evaluated at the 
projected distance p « 100 kpc of the M31-M33 field pre¬ 


dicts Nsni = 2.1 x 10 13 cm -2 , which implies an average 
total column density Nh = 6 x 10 17 (Z 0 /Z) cm -2 . Thus, 
if the CGM of M31 has a sub-solar metallicity, even if it 
is > 90% ionized, it would have a neutral component of 
similar magnitude to the average Nhi of the clouds over 
our field. 

The tight kinematic pattern of the M31-M33 clouds 
suggests that they are not spread along the several hun¬ 
dred kpc path through the M31 halo, so arguments about 
their origin in the M31 CGM are suggestive at best. The 


Lehner et al. (2015) measurements, however, allow the 
possiblity that the existence of clouds at p — 100 kpc 
from M31 does not require a major enhancement to the 
mass or density of its CGM, but could result from a re¬ 
structuring or phase change in material already present, 
triggered, perhaps by the pass age of a satellite as sug¬ 
gested by Shaya & Tully (2013), or by a concentration of 
mass in a plane of satellites. 

The so called “galactic fountain” model of gas accre¬ 
tion (e.g. Bregman [1980 ; Shapiro & Field 1976) states 
that supernovae kick material above the disk of a spiral 
galaxy, which then rains back down onto the disk as it 
cools. More recent work has shown that this fountain 
material can cause the surrounding halo material to con¬ 
dense and fall onto the disk as well (Fraternali & Binney 
2008| |Fraternali et al.|2015|). While the M31-1V133 clouds 
are too far away from either galaxy to be triggered by 
a galactic fountain, it is possible that hot gas could be 
triggered to condense by other external sources, such as 
the motion of a satellite galaxy through the CGM. Such 
a satellite galaxy would presumably need to contain at 
least some cold gas to serve as a trigger. It would be 
interesting to evaluate this mechanism in the context of 
the M31 CGM. 

7. SUMMARY 

Prompted by the discovery of extended regions of very 
faint HI that appear to form a partia l bridge connecting 
M31 and M33 (Braun & Thilker 2004), we have measured 
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21cm HI emission in a 12 square-degree region between 
M31 and M33 with the Green Bank Telescope at 9.'1 
angular resolution reaching 5cr limits on Nhi of 3.9 x 10 17 


, —2 


.-1 


also made at 18 locations on a grid to the north of M31. 

The new data confirm and extend the basic picture 
derived from the preliminary data presented in Paper I: 
in this region between M31 and M33 the HI is largely 
if not entirely contained in discrete neutral clouds that 
each have Mhi reaching a few 10 5 M©, lying at a pro¬ 
jected distance ~ 100 kpc from M31. We do not find 
any evidence for a more diffuse component of HI , and 
attribute our claim for this in Paper I to systematic in¬ 
strumental baseline effects at the level of a few mK. We 
measure only 63% of the HI mass found by BT04 in this 
region. While we present a possible explanation for this 
discrepancy (j |3.2 ), its origin is uncertain. The clouds ap¬ 
pear to be spatially and kinematically independent from 
each other and can have velocities that differ by > 100 
km s -1 over projected distances ~ 10 kpc. Our H I mass 
limits of ~ 10 4 M© are lower, by an order of magnitude 


or more, than other surveys of gala xy groups (e.g Auld 
et al.||2006 Chynoweth et al.||2008 ). Thus these objects 


may represent a new, previously undetected population. 

The clouds have a dynamical mass nearly a thousand 
times their HI mass and strong limits on any stellar 
component, making it unlikely that they are part of the 
dwarf galaxy system of the Local Group. Indeed, dwarf 
galaxies near large s pirals i n th e Local G roup comp letel y 
lack detectable HI ( Grcevich fc Putman||2QQ9 Spekkens 


|et al.||2Q14| |Westmeier et al.||2U15[ K. Beaton private 
communication). The clouds have kinematics more sim¬ 
ilar to the systematic velocity of the galaxies than to 
the HVC system of M31 and M33, but the clouds have 



if the M31-M33 direction is indeed a region of enhanced 
total mass. 

Our results to the north-west of M31 are still too in¬ 
complete to determine if we have detected anything like 
the population of clouds that exist between M31 and 
M33, but the data do suggest that the HVC population of 
M31 extends to p « 100 kpc in the north-west, much fur¬ 
ther than previously known. A similar extension to the 
south-east is not observed. The detection of HI outside 
the BT04 contours at S = +47° would imply that we are 
seeing very small angular sized features that are beam- 
diluted in the BT04 measurements. Complete mapping 
of this area with the GBT is underway to resolve the 
discrepancy. 


We thank Rachel Beaton for sharing results of her 
GBT survey of M31 dwarf galaxies before publication. 
S.A.W. acknowledges partial support from the student 
observing support grant (GSSP11-012) provided by the 
NRAO. D.J.P. and S.A.W. acknowledge partial support 
from NSF CAREER grant AST-1149491. 
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